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ABSTRACT 


The MIRI Si:As IBG detector arrays extend the heritage technology from the 
Spitzer IRAG arrays to a 1024 x 1024 pixel format. We provide a short discussion 
of the principles of operation, design, and performance of the individual MIRI 
detectors, in support of a description of their operation in arrays provided in an 


accompanying paper (Ressler et ah (2014)). We then describe modeling of their 
response. We hnd that electron diffusion is an important component of their 
performance, although it was omitted in previous models. Our new model will 
let us optimize the bias voltage while avoiding avalanche gain. It also predicts 
the fraction of the IR-active layer that is depleted (and thus contributes to the 
quantum efficiency) as signal is accumulated on the array amplifier. Another 
set of models accurately predicts the nonlinearity of the detector-amplifier unit 
and has guided determination of the corrections for nonlinearity. Finally, we 
discuss how diffraction at the interpixel gaps and total internal reflection can 
produce the extended cross-like artifacts around images with these arrays at short 
wavelengths, ~ 5 /rm. The modeling of the behavior of these devices is helping 
optimize how we operate them and also providing inputs to the development of 
the data pipeline. 


Subject headings: instrumentation: detectors; space vehicles: instruments 


^ Steward Observatory, 933 N. Cherry Ave, University of Arizona, Tucson, AZ 85721, USA 

^Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr. Pasadena, CA 91109, 
USA 

^Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 

^Laboratoire AIM Paris-Saclay, CEA-IRFU/SAp, CNRS, Universite Paris Diderot, F-91191 Gif-sur- 
Yvette, France 

^I. Physikalisches Institut, Universitat zu Kdln, Ziilpicher Str. 77, 50937 Koln, Germany 
®NASA Ames Research Center, M.S. 245-6, Moffett Field, CA 94035, USA 
^RALSpace, STFC, Rutherford Appleton Lab., Harwell, Oxford, Didcot 0X11 OQX, UK 






2 


Introduction 


The detectors of choice for the 5-28 /im range are arsenic-doped silicon impnrity band 
condnction (Si:As IBC) devices. They have extensive space flight heritage, for example, 


arrays of these devices were nsed in all three Spitzer instrnments (IRAC: Hora et al. (2004a); 


IRS: 

van Cleve et ah 

(1995 

), 

Houck et al. 

(2004 

); and MIPS: Gordon et ah 

(2004)), in 

WISE ( 

Mainzer et al. 

2008 

, in MSX ( 

Mill et al. 

1994 

) and in Akari ( 

Onaka et al. 2007). 


The focal planes on these missions have demonstrated high detective qnantnm efficiencies, 
low dark cnrrent and relative freedom from other spnrions signals, excellent photometric 
performance, and resistance to the effects of cosmic radiation. Similar detector arrays were 
selected for the Mid-Infrared Instrnment (MIRI) on JWST. They have been constrncted 
at Raytheon Vision Systems, specihcally based on the experience at this snpplier with the 


detector arrays for IRAC (Wn et al. 1997). The readont integrated circnit (ROIC), electronic 
data chain, mechanical arrangement, and operation of the MIRI arrays are described in 


Ressler et ah (2014), hereafter Paper VIII; while Gordon et al. (2014; Paper X) describe 


plans to calibrate the detector signals, inclnding removal of the artifacts discnssed below. 
The cnrrent paper discnsses the individnal detectors themselves.After an introdnction to 
the basic detector properties in Section 2, in Section 3 we describe work to nnderstand 
the detector performance from hrst principles as a way to help optimize the instrnment 
performance and to snpport the development of the data redaction software. 


2. Detector design and performance 
2.1. Detector design 

The structure of the detector layers (alternately referred to as Blocked Impurity Band— 
BIB—or Impurity Band Conductor—IBC device, depending on the vendor) is shown in Fig¬ 
ure 1. The pixels are manufactured on a wafer that provides a transparent silicon substrate, 
on which is produced a buried transparent contact and then an infrared-active detector layer 
of thickness 25 - 35/rm, relatively heavily doped (with arsenic for the MIRI devices) to pro¬ 
duce an impurity energy band. This layer is followed by an intrinsic blocking layer about 
3 - 4/im thick, with the second contact on it. This contact dehnes the pixel and is indium- 
bump-bonded to the input of the ROIC ampliher (see Figure]^. The bias voltage on the 
pixel is established between its individual contact and the buried common one whose bias 
is maintained through a V-shaped etched trough, aluminum-coated to make it conductive, 
and placed to one side of the array (Figure [^, although for illustration the hgure shows it 
in the middle. 
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Impurity band conductivity in the IR-active layer would produce unacceptable dark 
current if the band had electrical access to both detector contacts. However, the impurity 
band terminates at the intrinsic blocking layer, so leakage can only occur through thermal 
excitation up into the conduction band (this discussion is for an n-type dopant, e.g. arsenic). 
So long as the detectors are sufficiently cold (e.g., below 6.7K), this form of thermal current 
is very small. Because the intrinsic layer blocks the dark current, the doping in the infrared- 
active layer can be two orders of magnitude greater than in a bulk photoconductoij^ The 
highly doped, relatively thin IR-active layer overcomes the issues with bulk detectors of slow 
adjustment toward electrical equilibrium and of large cross section for ionizing radiation (for 


more discussion of these detectors, see, e.g., Petroff & Stapelbroek (1985); Rieke (2002)) . 


To collect the photo-generated charge carriers, the electrical conductivity of the IR- 
active layer must be suppressed to allow a field to be generated across it. When a positive bias 
is applied across the blocking layer to the buried contact, negative carriers are collected at 
the interface between the intrinsic blocking and IR-active layers and positive ones are driven 
away from this interface. Thus, a high-resistance region depleted of free charge carriers is 
produced in the IR-active material near the interface. However, beyond the depletion region 
there is virtually no field; this situation is shown schematically to the right in Figure 1 where 
the field is indicated by gray-shading and it dies away before reaching completely across 
the IR-active layer. Photoelectrons are produced throughout the IR-active layer; to gather 
them efficiently requires that it be fully depletecQ The minority, p-type impurities attach 
electrons and maintain a negative space charge in this region that tends to cancel the effect 
of the positive bias. The expression for the width of the depletion region shows this effect: 

1/2 


W = 


q Na I'^biasl + t 


B 


— ts 


( 1 ) 


where Na is the p-type impurity concentration, q is the electronic charge, Kq is the dielectric 
constant of silicon, Sq is the permittivity of free space, \Vuas\ is the bias voltage across the 


detector contacts, and ts is the thickness of the blocking layer (Rieke 2002). Equation (1) 
implies that large values of Na can be overcome with large bias voltages. However, if the 
voltage across the detector is made too large, the field near the blocking layer can become 
strong enough that the photoelectrons gain sufficient energy to free additional electrons, 
creating avalanche gain. This gain is undesirable because it is generated by a stochastic 
process and adds noise to the signal. Thus, there is an imperative in keeping the minority 


^That is, a detector consisting of an undifferentiated layer of photoconductive material placed between 
two electrical contacts. 

^Diffusion lengths of a few also allow some charge collection from regions adjacent to those that are 
depleted. 
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impurity concentration low to make w large at a modest level of Vuas- Complete depletion 
in the MIRI arrays without signihcant gain requires that the level of the minority impurities 
be kept below about 2 X 10^^ cm“^. That is, the impurity atoms must be at or below 1 in 
10^° silicon atoms. 


In Figure 1, light enters from the top and passes through the thick, inactive substrate 
layer of pure (i.e., high resistivity) silicon. After also passing through the transparent elec¬ 
trical contact, the photon is absorbed in the active layer of arsenic-doped silicon, as shown 
to the right in Figure 1. The arsenic concentration can be made high enough for almost 
complete absorption in the 10 - 25 /xm range, where the dominant losses occur at the trans¬ 
parent contact. However, the absorption cross section is wavelength dependent, and short 
wavelength photons may pass through the active layer, reflect off the metallized electrical 
contact, and be absorbed as they travel back out the active layer (or they may escape from 
the initial pixel, see discussion in Section 3.3.2 below). There are no physical boundaries 
between pixels in an IBC detector array — the detectors are all on a single slab of silicon. 
Pixels are dehned by the very high electric helds induced between the transparent contact 
separating the active and inactive layers and the contact at the indium bumps. The helds 
are so high that a created electron-hole pair is quickly swept out the nearest contact, so that 
interpixel diffusion is generally not an issue. 


The detector material for MIRI was grown at Advanced Technology Material, Inc. 
(ATMI), where precautions were taken to achieve high purity levels. Spreading resistance 
measurements suggest that the minority (acceptor) impurities were held to a density of 2 X 
10^^ cm“^ or less (Love et al. 2005). Two designs were implemented. The “baseline” arrays 
have an arsenic doping level of 7 X 10^^ cm“^ and an infrared-active layer 35 pm thick, 
while the ’’contingency” arrays are similar but with the IR-active region reduced to 30 pm 
in thickness and an arsenic concentration of 5 X 10^^ cm“^ . Both types have blocking layers 
of 4 pm thickness. The contingency material was designed to be virtually identical to that 
used in the IRAC arrays and thus to have a high probability of reproducing their perfor¬ 
mance. The baseline design was developed for the potential of higher quantum efficiency 
where the absorption cross sections are small (e.g., at wavelengths short of 8 pm) using a 
thicker infrared-active layer. The price is a higher dark current: about 0.2 e/s at the 6.7K 
operating temperature, about three times higher than for the contingency material. 


To manufacture detector arrays, the detector wafers were diced and patterned with 
indium bumps, allowing them to be hybridized onto the readouts (Figure]^. Each detector 
is connected to the input of a source follower ampliher as shown in Figure The array 
format is 1024 X 1024 pixels with a pitch of 25 pm. The hybridized arrays are then anti¬ 
reflection (AR) coated with either of two single layer AR thicknesses, one optimized for 6 
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/im and the other for 16 fim. These steps complete the construction of a MIRI Sensor Chip 


Assembly (SCA). The properties of the MIRI Si:As IBC SCAs are described by Love et ah 


(2004 

); Love et ah ( 

2005) 

Love et ah 

(2006 

) and 

Ressler et ah 

(2008 

), and have a heritage 


in both the detectors and ROIC to the arrays developed for IRAC (Wu et ah 1997; Hora 


et ah 2004a). 


2.2. Performance 


We introduce a few basic properties of the MIRI detectors to motivate the discussion of 


the operation of the devices in Section 3. The basic parameters are summarized in Ressler 


et ah (2008) and in Table 1. 


2.2.1. Radiometric properties 

Figure shows the response of the baseline detector material as a function of bias 
voltage at a wavelength near 13/rm. The general characteristics are a rise in response from 
zero bias to a bias of ~ 1.2 V as the increasing voltage depletes more of the infrared-active 
layer, followed by a leveling off for biases of 1.2 V < V < 2.5 V , with finally the beginning 
of an upturn in the response above 2.5 V as avalanche gain begins to become signihcant. As 
the detector photocurrent is collected on the input of the unit cell amplifier (Figure 3), the 
net voltage across the detector decreases and the response to further illumination is reduced 
according to the relation illustrated in Figure 4. As a result, the integration ramps have a 
smooth departure from linearity; just before the signals saturate, the reduction is about 10% 
relative to a linear extrapolation of the small-signal behavior. 

Figure]^ shows the measured responsive quantum efficency (using a large process control 
detector) as a function of wavelength, and the values computed from the same measurements 
assuming single layer anti-reflection coatings optimized for 6 and 16/im. 


2.2.2. Imaging properties 

The response of the detectors in an array is uniform, with pixel-to-pixel variations of 
order 3% rms. The best arrays have a small proportion of inoperative pixels (either hot or 
dead), of order 0.1%. 

The detectors have a complex latent image behavior. After removing a bright source. 
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there is an initial fast decay with a time constant of abont 8 seconds. After this signal has 
died away, there are latent images at the initial level of abont 1% of the initial image, in 
observations with a constant 3-second readont cadence. The fading of these images shows 
mnltiple time constants, as is shown in Fignre Extrapolating the three exponential hts 
in the hgnre back to the time the sonrce was removed indicates an initial latent image of ~ 
1% of the illnmination, in agreement with measnrements. Therefore, the three-exponential 
£t is not missing any significant faster effect. However, it appears that the details of the £t 
change according to the strength of the initial signal; it is not adeqnate to simply scale the 
response from one signal strength to another to £t the latent image behavior. 


The primary effect of the latent images is a depression of the reset level that is slowly 
restored to its pre-exposnre valne. The resnlt is that the hts to the integration ramps 
are steepened, creating an apparent positive image that decays in strength as the depression 
fades. The latent image behavior is similar to that observed previonsly with similar detectors. 
For example, the primary latent decay time constant for the MIRI array of 8 seconds is very 
similar to the 12 ± 5 seconds reported for the IRS/MIPS arrays (Gordon et ah 2005). The 
MIPS arrays also had mnltiple time constant latent decay patterns, snch as long-life latent 
images after hard satnration (Gordon et ah 2004). The IRAG arrays had similar behavior 


(Hora et ah 2004a). 


The imaging properties of MIRI were measnred mnltiple times dnring the bnildnp of 
the instrnment modnles and then in the test of the hight model prior to delivery, and hnally 
in the Integrated Science Instrnment Modnle (ISIM) test post delivery. The most precise of 
these measnrements in terms of the array performance were condncted at the Gommissariat 
a I’energie atomiqne (GEA) with jnst the imager, and illnminated by a sonrce ontside the 


cryostat with extensive filtering to control the backgronnd emission (e.g., Ronayette et al. 


(2010)). The point spread fnnction (PSF) measnrements ntilized a micro-stepping strategy 


so that many positions of the sonrce were recorded, on centers smaller than the pixel pitch of 
the MIRI array. These measnrements were then converted to a high-resolntion PSF image. 
The images at the longer wavelengths are as expected (Fignre [^, showing excellent imaging 
properties from the array. However, at 5.6 /im there is an additional cross-like imaging 
artifact as shown in Fignre Similar artifacts were seen with the IRAG arrays (Fignre [^. 

At all wavelengths, there is a low level of pixel-to-pixel crosstalk, which we have mea- 


Finger et al. 

(2(1(15 

), 

Regan & Bergeron 

(2012 

), 

Rieke & 


Morrison (2012)), inclnding antocorrelation, analysis of cosmic ray hits, and of hot pixels. 
All measnrements indicate a level close to 3% for the crosstalk to the fonr adjacent pixels 
aronnd the one receiving signal. The dominant canse of this behavior is interpixel capac¬ 
itance, althongh there may be secondary contribntions from electron diffnsion and optical 
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effects (Regan & Bergeron 2012 Rieke & Morrison 2012). 


3. Modeling 


3.1. Responsivity 

The theory of silicon impnrity band condnction detectors (Si:X IBC) was developed in 


a series of papers pnblished in the 1980s (Petroff & Stapelbroek (1984, 1985); Szmulowicz & 


Madarsz (1987); Szmulowicz et ah (1988)). However, those studies were of detectors with 
minority impurity concentrations of 10^^ - 10^® cm“^, majority impurity concentrations of 
about 5 X 10^^ cm“^, and infrared-active layers 15 to 25 /im thick. The MIRl detectors differ 
signihcantly with regard to the hrst and third parameters. We have therefore generated a 
new set of detector models, accounting for performance aspects that were not apparent in 
the older analyses. 

The model for detector responsivity will be illustrated for the baseline detectors. The 
basic input data are contained in Figure]^ and the parameters of the detectors are described 
in Love et ah (2004) and Love et ah (2005). The detailed theory of responsivity vs. bias of 


IBC detectors was discussed in Petroff and Stapelbroek (1984, 1985). In our model, we used 


equation (39) in Szmulowicz & Madarsz (1987) to quantify the avalanche gain as a function 


of depth in the detector. The gain is largely controlled by the characteristic held strength, 
Ec- The adopted value has a strong inhuence on the upturn in the responsivity curve for 
biases above 2.5V in Figure]^ we leave it as a free parameter in our model to allow accurate 
htting of this behavior. 

Previous models have assumed that the signal-producing region of the detector includes 
only that part of the IR-active layer that is depleted of charge carriers by the applied electric 
held as in equation (1) (e.g.,Petroh & Stapelbroek (1985); [Szmulowicz fc Madarsz (1987)). 
Photoelectrons produced in the undepleted (neutral hereafter) region are assumed to recom¬ 
bine without contributing to the signal. This assumption is difficult to reconcile with the 
behavior of the MIRI detectors. At the wavelengths near peak absorption, the absorption 
cross section is of order 1.3 X 10“^® cm^. Because of the relatively large thickness of the 
IR-active layer and the high concentration of As, the absorption length is only about 10 
fim. In our model, at 0.5V bias the neutral layer is about 18 fim thick and for our back- 
illuminated detectors the absorption occurs at the side of the IR-active layer opposite the 
depleted region, leaving a gap of about 8 /xm between the depleted region and where the 
photons are absorbed. Thus, one would conclude that there would be very little signal at 
biases around 0.5V, in disagreement with the slow reduction in response with reduced bias 






































in Figure 


The missing ingredient is diffusion of the free charge carriers from their places of origin. 
There is no analytic form for the general case of diffusion of charges within a plane-parallel 
detector, so we have used the expression for the case where the absorption occurs all at the 


entrance to an infinite planar detector and diffuses across a region of thickness c (Holloway 


1986): 


S 


2b 


( 2 ) 


where L is the diffusion length and b is an efficiency term to allow for the input flux to the 
diffusion region. We have applied this model iteratively across the neutral region in steps of 
1 micron to simulate the situation in the MIRI detectors, where the absorption can extend 
well into this region. 


Figure shows the results from a model for the baseline array where the impurity 
concentration, characteristic field strength, and diffusion length were varied to match the 
voltage dependence of the response. The parameters for this model are Ec = 8700 V/cm, 
cm“ 

-2 


Na = 1.5 X 10^^ cm and L = 2.5 pm. The absorption cross section was set to 5 X 


cm ^ (Geist 1989). The derived minority impurity concentration is consistent with 


10-16 

the spreading resistance measurement (Love et ah 2005). The corresponding value for the 
recombination time is about 2 X 10“^ s, assuming a mobility of 3 X 10^ cm^ s“^ from 
the I. F. Ioffe compilation (Ioffe 2014). The model indicates that the avalanche gain at 
the selected bias voltage of 2.2V is only a few %. A second model was produced for the 
contingency array data, as shown in Figure 10 The parameters are Ec = 8900 V/cm, 
Na = 2.3 X 10^^ cm“^, and L = 2 pm. Except for the minority impurity concentration, 
they are in good agreement with the parameters for the baseline array, which supports the 
validity of the model. Unfortunately, it appears likely that the minority impurities are not 
as well controlled for this material. The avalanche gain at a bias of 2.2V is indicated to be 
about 10%. 


Although both models provide reasonably good fits to the observed trend of responsivity 
vs. bias voltage, neither is perfect. An obvious cause of discrepancies is that the models 
assume a single value of the minority impurity concentration throughout the IR-active layer. 
In practice, it is likely that the impurity concentration has gradients. Even more fundamen¬ 
tally, a 10 pm thick layer under a pixel will have only about 10,000 impurity atoms; that 
is, the rms statistical fluctuations will be at the 1% level, sufficiently large to influence the 
responsivity behavior. The changes in the responsivity of a given pixel due to just the statis¬ 
tical fluctuations in the minority impurity concentration under that pixel are large enough 
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to contribute significantly to the rms pixel-to-pixel responsitivity variations. 


3.2. Nonlinearity 


Although nonlinearity in near-infrared photodiode arrays is dominated by changes in 
the pixel capacitance, the small fraction of the node capacitance due to the detectors in the 
MIRI arrays demonstrates that this mechanism is unimportant for them. To be specific, the 
capacitance at the input to the unit cell buffer amplifier is 28.5 fF (McMurtry et ah 2005), 


the detector is a simple plane-parallel 1.1 fF capacitor, and the bump bonds may contribute 
an additional ~ 4 fF to the node capacitance ( Moore [2005 ). The nominal capacitance at the 
integrating node is therefore 33.6 fF, of which only about 3% is contributed by the detector. 


The array amplifiers are linear to within about 1%, as measured by injecting signals into 
the integrating node to calibrate the gain as a function of signal strength. The dominant 
cause of nonlinearity is the reduction in detector responsivity as signal is accumulated on 
the amplifier integrating node. When the detector bias is reset, the voltage across it will 
be set accurately to Vbias (see Figure 3) because its effective resistance of ~ 5 x lO^^hl 
(estimated from dark current measurements) is far larger than the “on” resistance of the 
reset FET. Once the reset FET is turned off, as charge accumulates on the amplifer input 
node capacitance it reduces the net bias voltage on the detector. The result is that the 
responsivity declines as shown in Figure 

A model based on this behavior is able to reproduce the array nonlinearity accurately. 
The performance is illustrated in Figure which shows the ratio of a model signal assuming 
perfect linear response to the measured signal corrected for nonlinearity. In computing this 
ratio, the only free parameters are the charge accumulation rate (the slope of the integration 
ramps) and the gain of the detector/amplifier system (that is the DN of output per volt 
change at the integrating node, which determines the relation of the DN of the signal to the 
position on the responsivity curve in Figure 4). As shown in Figure 11, the best fit to the 
shape of the nonlinearity (that is where the residuals relative to a ratio of 1 are smallest) is 
for a gain of about 35000 DN/V. The voltage gain of the system has been measured directly 
to be 38300 DN/V, determined by injecting a signal into integrating node and measuring the 
output (Paper VIII). That is, the agreement between the value deduced indirectly through 
the model for nonlinearity and the measured one is excellent. Combined with the node 
capacitance, discussed above, the pixel gain is 5.5 e/DN, with an uncertainty of about 10%. 


As charge is collected, eventually the detector bias will fall low enough that the depletion 
depth of the IR-active layer will be reduced. According to the detector model described in 
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Section 3.1, the baseline array is no longer fully depleted at biases of 1.5V and lower, that is 
after collection of a total signal of about 27,000 DN or 150,000 electrons. However, the loss 
of depleted depth is fairly slow; at IV, about 75% of the IR-active region is still depleted. 
A bias of one volt corresponds to an integrated signal of 46,000 DN or 250,000 electrons; 
it is a useful limit for ’’full well” and in fact the MIRI warm electronics saturates not far 
above it. Nonetheless, there are observable consequences from the loss of full depletion: 1.) 
for large signals, the detective quantum efficiency of the detectors will be reduced; and 2.) 
because of the different absorption prohles in the IR-active layer for different wavelengths, 
the nonlinearity will have a wavelength dependence. 


3.3. Imaging properties 

3.3.1. Latent images 


Unfortunately, there appears to be little previous analysis of potential causes of latent 
images in these arrays, so we will make some suggestions without coming to firm conclu¬ 
sions. Low levels of trapped charge could result from impurities in the blocking layer, which 
might cause weak latent images at all levels of illumination. Another possibility is dielectric 
relaxation - the slow adjustment of the electrical equilibrium of some portion of the detector, 
basically because of its long RC time constant. Stapelbroek (private communication) has 
calculated that the time constant for this process in the undepleted region of the detectors 
is roughly appropriate to lead to the more rapidly decaying latent images. [Smith et ah 


(2008) propose that latent images in photodiodes can arise because the depleted region at 


the junction shrinks as the detector is debiased due to accumulation of signal, allowing trap¬ 
ping of photoelectrons. Once the bias is re-established, the traps lie within the depletion 
region; as the trapped electrons are slowly released, they appear as a signal. Although this 
argument cannot be applied directly to Si;As IBC detectors, a similar situation may exist 
if sufficient signal is gathered to produce a significant undepleted region (from the above 
section, a signal of order 200,000 electrons would suffice). A number of these processes could 
operate together. Thus, why the latent image behavior might be complex is not difficult 
to understand, but a systematic set of measurements as a function of detector temperature, 
photon flux, and perhaps wavelength is needed to obtain sufficient information to constrain 
the possibilities. 




3.3.2. Image ’halos’ at short wavelengths 


Observations with IRAC (in spectral bands at 5.8 and 8 /rm) showed that a signihcant 
fraction of light was being detected in what was described as a halo around the core image 
(Fazio et ah 2004; Hora et ah 2004b; Pipher et ah 2004; Figure 8). This behavior was 
also seen in laboratory measurements (C. McMurtry, private communication), where the 
phenomenon appears to be more accurately described as light scattered along the rows and 
columns of the array (e.g., a cross geometry centered on the image core). The behavior is 
specific to operation of this type of detector array at wavelengths between 5 and 10 /xm (and 
presumably at shorter ones, but the devices are seldom used short of 5 jam). The clearest 
demonstration that this behavior is absent at longer wavelengths is the excellent agreement 


between the point and extended source calibrations at 24 /xm for MIPS (Cohen 2009). 


Figures 7 and 8 demonstrate that the MIRI arrays have very similar imaging performance to 
those in IRAC and MIPS, i.e., clean images at wavelengths longer than 10 /xm, but cross-like 
artifacts at shorter wavelengths. 

An explanation of this behavior is based on the geometry illustrated in Figure 2. The 
detectors are fabricated on a silicon substrate ~ 500 /xm thick that is transparent in the mid- 
IR. Light enters this substrate at the top of the figure. Even with AR coatings, reflectivities 
in the 5 - 8/xm range are of order 8% (Hora et ah 2004b Love et al. 2005). The mid- 


IR photons are absorbed in the IR-active layer. Between it and the bulk of the substrate 
there is a thin transparent buried contact, which can also reflect some of the light; nominal 


reflectivities can be up to 20% (Petroff & Stapelbroek 1985), or 25% (Szmulowicz et al 


1988). The reflection losses are also apparent from Figure 5, since the responsive quantum 


efficiency never gets above about 80% even at wavelengths where the absorption in the 
IR-active layer should be complete and the AR coating highly efficient, e.g., ~ 16 /xm. 

For the 7 x 10“^^ cm“^ arsenic concentration in the IR-active layer of the MIRI baseline 


detectors, the imaginary component of the index of refraction is (Woods et al. 2011): 


k = 5.69 X 10 


-3 


A 


7 /xm 


( 3 ) 


Here A is the wavelength. The corresponding linear absorption coefficient is: 

“ = (pL) 

(This formulation is not applicable at the longer wavelengths, e.g., >15 /xm.). For the 35 
/xm active layer thickness of the MIRI baseline detectors, the absorption at normal incidence 
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is 17%, 65%, and 95% at wavelengths of 5, 12, and 20 /xm for a single pass throngh the 
IR-active layer. Most of the light not absorbed will be reflected by the pixel contact to make 
a second pass, so the total absorption will be abont 28%, 83%, and 100% respectively at 
these three wavelengths. At short wavelengths, a signihcant amonnt of light snrvives after 
enconntering the active layer. To understand the cross artifact, we need to understand the 
fate of this light. 


Simple arguments show that the cross artifacts must be associated with the interpixel 
gaps (see Figure]^. Optically, these gaps will act as narrow slits, with widths comparable to 
the wavelength of the light. Therefore, we cannot use a geometric approach to understand 
their effects, since they will spread the light by diffraction. Traditionally, diffraction is 
treated in terms of the effects on the light transmitted through a slit. However, Babinet’s 


Principle (Hecht (2001), Section 10.3.11) states that the Fraunhofer diffraction patterns 


of complementary masks are identical, that is the equivalent diffraction pattern will be 
obtained if the slit is replaced by an opaque occulter with free transmission in the surround. 
This situation obtains for the reflected light around an interpixel gap in the detector arrays 
(with the addition of a fold of the beam through the reflection). We can estimate the 


angular dependence of the diffracted light in terms of a mm-wave analog from Lee et ah 


(2009 who illuminated slits 10 /rm thick with light of wavelength between 6000 and 150 


/im, comparable to the ratio of wavelength to slit thickness of about 50 for the contacts on a 


Si:As IBC array. Lee et ah (2009) hnd transmissions approximating the geometric thin slit 


case for frequencies up to a cutoff of c/2a where c is the speed of light and a is the slit width. 
In the case of the Si;As IBC arrays, the physical interpixel gaps (slit widths) are effectively 
magnihed by the refractive index of silicon, so they can be considered to be 27 pm (IRAC) 
and 10 pm (MIRI), and fall into this pseudo-geometric regime. Thus, far from the gaps the 
illumination pattern should approximate classical Fraunhofer diffraction. 


Therefore, a single-slit diffraction pattern will be imposed on the reflected beam equiv¬ 
alent to that of a slit with the geometry of the interpixel gap. Some of the diffracted light 
will be reflected off the buried contact and the input surface of the detector to be detected 
removed from the pixel where it entered the array. The spreading will be perpendicular to 
the slits, which explains why the artifacts are along rows and columns of the array. However, 
simple diffraction and reflection do not account for the extent of the cross artifact. A more 
detailed model shows that some of the light is diffracted beyond the angle for trapping by 
total internal reflection in the detector. In addition, some of the light escapes the detector 


^ Other studies of the effects of narrow slits on optical wavelength light tend to use slits with thick¬ 
nesses larger than the photon wavelength and the emergent light angular distribution is strongly affected by 
waveguide effects. 
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wafer through the interpixel gaps into the space between the detector and the readout circuit 
(see Figure 2). This region acts like a lossy integrating cavity, and some of the light can 
re-enter the detector wafer through other interpixel gaps. Again, this light may be diffracted 
to angles where the photons are trapped by total internal reflection. Because the entire 
wafer is inefficient at absorbing photons in the 5 /im region, those that are trapped can 
travel large distances before they are absorbed and produce signals. This behavior can ac¬ 
count qualitatively for the cross image artifacts; a more detailed discussion will be published 
elsewhere. 


4. Conclusion 

The Si:As IBC detector arrays for MIRI represent the next step in the development of 
detectors for the 5-28 /rm spectral range, building on the heritage from the detectors used 
in the three Spitzer instruments, WISE, MSX, and Akari. They have a larger 1024 X 1024 
pixel format, and show good quantum efficiency and low dark current. We have modeled 
the detector behavior to understand the avalanche gain characteristics, the behavior of their 
nonlinear response, the effect of large signals on their absorption efficiency, and the extended 
images found near 5 /rm. This work will help optimize the operating parameters for the 
detectors and provides inputs to the data pipeline development for MIRI. 
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Table 1. 


Detector Performance Parameters 


Parameter 

baseline array 

contingency array 

format 

1024 X 1024 

1024 X 1024 

pixel size 

25 /im 

25 /im 

IR-active layer thickness 

35 /im 

30 /im 

IR layer As doping 

7 X 10^^ cm“^ 

5 X 10^^ cm“^ 

read noise* 

14 e" 

14 e" 

dark current 

0.2 e“/s 

0.07 e-/s 

quantum efficiency** 

> 60% 

> 50% 

nominal detector bias*** 

2.2V 

2.2V 

well capacity 

~ 250,000 e- 

~ 250,000 e- 


* Fowler-eight sampling, used for comparison purposes; the read¬ 
out is normally operated in a sample-up-the-ramp mode. 

**At peak wavelength 

***Consisting of 2 V applied directly plus ~ 0.2 V from clocking 
signal feedthrough 
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Buried 

Transparent. 

Implant 

High Resistivity (Transparent) Substrate 

( 

r 

Detector Implant^ 

Active Detector //V Buried ^ 

^Layer (As-doped)/ . 




Oxide 
Indium Bump 


Al Contact to Buried 
Transparent Implant 


Detector Contact 


Fig. 1.— Basic operation of a Si:As IBC detector array. The left side shows the physical 
arrangement of the detectors. They are manufactured on a lightly doped (high resistivity) 
transparent substrate. A 25 - 35 /rm thick layer is relatively heavily doped with arsenic and 
absorbs the incident photons. An electric held is maintained across this layer by putting a 
voltage between the detector contact and the buried transparent common contact; connection 
to the buried contact is made through the V-etch (to one side of the detector array). The 
held is shown schematically by the gray-scale on the right pixel. An intrinsic layer blocks 
thermally generated free charge carriers in the arsenic impurity level from escaping from 
the IR-active layer and reaching the detector contact (because there is no corresponding 
energy level in intrinsic material), but allows the photoelectrons passage since they have 
been elevated into the conduction band. A photon is shown generating a photoelectron just 
at the edge of the depletion region where there is sufficient electric held to draw the free 
electron to the appropriate detector contact. When the array is hybridized to the readout 
wafer, contact to an output ampliher is made through an individual indium bump for each 
pixel. Figure based on that in Love et al. (2005). 
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Fig. 2.— Illustration of the geometry of an arsenic-doped silicon detector array. The light 
travels through 500 /im of pure silicon (which may be anti-reflection coated) and then a 
transparent contact into a thin active layer (25 pm for IRAC, 35 pm for MIRI baseline 
material) of arsenic-doped silicon, followed by an even thinner (4 pm for MIRI) blocking 
layer of pure silicon. Indium bump bonds, one per pixel, join evaporated contacts on the 
detectors and readout amplihers to convey the photo-electrons generated in the active layer 
to the input of the amplifiers. 
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Fig. 3.— A representative Si;As IBC detector array unit cell. The circuit is a simple 
source follower, where current from the detector (drawn as a diode to indicate its electrical 
assymetry) delivers charge to the gate of the analog MOSFET ampliher. Vout is proportional 
to the charge integrated onto the gate. The unit cell is addressed by Ti (column) and T 2 
(row), and T 3 is used to set the detector voltage back to the initial bias level. The design is 
based on that of the unit cell for the IRAC arrays. 
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Fig. 4.— Relative responsivity vs. bias voltage for the baseline material (Ressler et al. 


2008). Good control of the minority impurity in the IR-active layer is indicated both by the 


knee near 1.2V which indicates that full depletion of that layer is being approached at that 
voltage, and the slow onset of avalanche gain above 2.5V. 
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Fig. 5.— Measured responsive quantum efficiency of bare detector material (solid line). The 
dashed line is a computed result assuming the array has an antireflection coating applied 
optimized for 6 /rm, and the dotted line is for an AR coating optimized for 16 pm. This 


hgure is based on one in Ressler et ah (2008). 
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Fig. 6.— Recovery from a saturating signal. The signal was about 20,000 e/s and overfilled 
the ampliher wells by about a factor of two (i.e., about 500,000 total electrons); there is about 
a 3 minute interval between removing the signal and the start of data collection, to allow 
placing the detector in the dark. The signal (points) has been htted with three exponentials: 
dot-dash, timeconstant of 130 seconds; long dashed, time constant of 14 minutes; short 
dashed, time constant of 130 minutes. The two longer time constants are very uncertain 
because they depend strongly on the placement of the zero signal baseline. The solid cnrve 
shows the snm of the three fits. 
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Fig. 7.— Measured point spread function with a simulated JWST pupil, from left to right, 


for 7.7 /im, 12.8 /rm, and 15 /rm (Ronayette et ah 2010) 
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Fig. 8.— Left: measured point spread function for the MIRI imager at 5.6 pm (Ronayette 


et al. 2010). Right: image with IRAC at 5.8pm of the bright protostar S140 (Pipher et ah 


et al. 

2010 

2004 

). The 


The cross response is an artifact resulting from the low absorption cross section for 


these photons in the IR-active layer. 
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Fig. 9.— Model responsivity vs. bias (line) compared with measurements (points) for the 
baseline detector material. 



relative response 
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Fig. 10.— Model responsivity vs. bias (line) compared with measurements (points) for the 
contingency detector material. 
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DN 


Fig. 11.— Departure from linearity for different assumptions about the system gain. A zero 
point of 6000 DN has been subtracted from the integration ramp before taking the ratio 
relative to the predicted purely linear behavior. 






